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Abstract. Noting that the number of gluons in the hadron wave function is discrete, and their 
formation in the chain of small x evolution occurs over discrete rapidity intervals of Ay ~ l/(X s , 
we formulate the discrete version of the Balitsky-Kovchegov evolution equation and show that its 
solution behaves chaotically in the phenomenologically interesting kinematic region. 
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EVOLUTION AS A DISCRETE PROCESS 

The color field of an ultra-relativistic hadron is a quasi-classical non-Abelian 
Weizsacker- Williams field yj, yfl. It emerges when the occupation number of the 
bremsstrahlung gluons emitted at a given impact parameter exceeds unity and eventu- 
ally saturates at ~ 1/ a s . It has been argued that in a big nucleus, such that a s A 1 ^ 1, 
and not very high energies the mean-field treatment is a reasonable approximation to 
the evolution equations. 

As a result of broken scale symmetry of QCD there exist the dimensional scale A 
which is the infrared cutoff on the gluon's momenta. An introduction of an infrared 
cutoff A on the momentum of the emitted gluons amounts to imposing the boundary 
condition. This is equivalent to the quantization of the gluon modes in a box of size L ~ 
A 1 , in which case the spectrum of the emitted gluons and their number become discrete. 
The formation of a gluon occurs over a rapidity interval of Ay ~ l/ot s . Therefore, 
the evolution in rapidity can be considered as a discrete quantum process, where each 
subsequent step occurs when Ay . ~ 1 01]. 

Assuming that Ay a s is a certain number for all steps in evolution process neglects 
a stochastic nature of quantum evolution. Full treatment of the discrete BK equation 
requires taking these effects into account. However, unfortunately BK equation is known 
to resist all attempts of analytical solution, and our hope at present is to develop a 
meaningful approximation. Thus, we suggest an approximation in which the gluons are 
emitted over a fixed "time" defined by a s Ay = C with C = 1 . To justify this assumption, 
let us note that BFKL takes into account only fast gluons, i.e. those with C ~ 1. It is 
beyond the leading logarithmic (LL) approximation to take into account slow gluons. 
Moreover, it is known that an account of NLL corrections effectively leads to imposing 
a rapidity veto ||4|] on the emission of gluons with close rapidities, which restricts 




production of gluons with small C (this is due to an effective repulsion between the 
emitted gluons induced at the NLL level). Therefore C is bounded from below by a 
number close to one. On the other hand the probability that no gluon is emitted when C 
becomes larger than one is very small if we choose the high density initial condition, 
such as the one given by the McLerran-Venugopalan model[l]. Therefore, C takes 
random values around 1, but the effective dispersion can be expected quite small. 



DISCRETE EVOLUTION EQUATION 

Equation which describes the gluon evolution in the high parton density regime of QCD 
is the Balitsky-Kovchegov equation This equation is formulated for the forward 

scattering amplitude of a color dipole of transverse momentum k and rapidity y to scatter 
of a big nucleus at impact parameter b. Assumption that the dispersion of the dipole 
transverse momenta with y is a slow process one can develop a diffusion approximation 
to the BK equation. This approximation is meaningful in the saturation region. The 
discrete version of BK equation takes the form (see also 0]) 

N n +l(k,y) = (l+x(Yo))N n (k,y) - N 2 n (k,y) , (1) 

where we introduced the discrete variable n to enumerate emitted gluons and x(Yo) = 
41n2 = ft) — 1. It is convenient to re-scale the scattering amplitude § n = x(Yo)N„ so that 
the corresponding continuous amplitude is normalized to Nnik.y) < 1: 

n+ i = (0(j) n - (co - 1)0„ 2 . (2) 

Numerical solution to the discrete BK equation is shown in Figures [l]-|4| 

Let us now consider how does the evolution proceeds for various values of CO. In Fig.[l] 
the case of 1 < CO < 3 is shown. In that case, for any k there is one stable fixed point at 
§ n = 1 and one unstable fixed point at § n = 0. The fixed points are determined from the 
condition (j) n+ i = ty n . 

In Fig. El we consider the case 3 < CO < 3.442 The point (j) n = 1 ceases to be 

unstable. Instead two new stable points appear. These can be determined from the 
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FIGURE 2. Discrete BK equation at co = 3.1. Different lines from left to right correspond to n = 
1,4,7, 10 (black, red, green, blue). 



FIGURE 3. Discrete BK equation at CO = 3.495. Different lines from left to right correspond to 
n = 1,4,7, 10 (black, red, green, blue). 



condition § n +2 = § n - The described multiplication of stable points is called in general 
bifurcation. In our particular case it is referred to as the period doubling scenario. 

When 3.442... < co < 3.56... there are four new fixed points, see Fig. It is 
important to emphasize that at n — > °° the value of (j) n settles to a given set of fixed points 
(specified by the value of co) independently of the initial condition. In other words, 
at very high energies the scattering amplitude in the saturation regime kj < Q s (n) is 
independent of the initial condition. 



The period doubling proceeds at increasingly smaller increments of CO until the accu- 
mulation point C0f = 3.569. . . known also as the Feingenbaum's number. At this point 
there is no more universal limiting behavior at large n. Instead small change in the initial 
condition leads to large change in the final state. The onset of this chaotic behavior can 
be observed in Fig. HJ Note now that at first, cObfkl = 1 +41n2 = 3.77 > COp and at 
second, CObfkl is the absolute minimum of the function 1 +x(Y)- Thus, we conclude 








ONSET OF CHAOS 
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FIGURE 4. Discrete BK equation at 03 = 3.11. Different lines from left to right correspond to n = 
1,4,7, 10 (black, red, green, blue). 



that the evolution of the scattering amplitude at high energies in the saturation region 
may be chaotic. At the same time, it is evident from Figs. [TJ-|4| that the high kj tail of 
the scattering amplitude which describes the perturbative regime is not affected by the 
peculiar behavior of the discrete equation. 

By averaging over all events one can define the mean value of the scattering amplitude. 
However, this procedure hides a lot of interesting physics. The most obvious example 
of this is diffraction, which measures the strength of fluctuations in the inelastic cross 
section. Figs. [TJ-|4| imply that diffraction is a significant part of the total inelastic cross 
section at very high energies, and is universal (independent of the properties of the 
target). 

The model used in this letter is admittedly oversimplified: we neglected the diffusion 
in transverse momentum, stochasticity of gluon emission and the dynamical fluctuations 
beyond the mean field approximation. Nevertheless, we hope that at least some of the 
features of discrete quantum evolution at small x will survive a more realistic treatment. 
The chaotic features of small x evolution open a new intriguing prospective on the studies 
of hadron and nuclear interactions at high energies. 
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